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. - . We apply the quantum Hamilton-Jacobi formalism, naturally defined in the com- 

, plex domain, to a number of complex Hamiltonians, characterized by discrete parity 
O ■ 

psj ■ and time reversal (FT) symmetries and obtain their eigenvalues and eigenfunctions. 

' Examples of both quasi-exactly and exactly solvable potentials are analyzed and the 

^ ! subtle differences, in the singularity structures of their quantum momentum func- 
tions, are pointed out. The role of the PT symmetry in the complex domain is also 

y—^ ' illustrated. 

> : 

vn ■ 

O ■ I- INTRODUCTION 

m : 
o . 

^ I Complex Hamiltonians possessing real eigenvalues have attrac ied considerable attention 

' in the current literature S H H H H Q • These quantal systems, 

apart from being counter intuitive, are not well understood because of their recent origin. 
■ These Hamiltonians are characterized by discrete parity and time reversal symmetries. In 

a ■ 

^ ' the case when the wave functions are also PT symmetric, the eigenvalues are real and the 



violation of PT symmetry by the wave function leads to eigenvalues, which occur in complex 
conjugate pairs. Besides identifying new Hamiltonians belonging to this class, the role of 



' various discrete symmetries is also under current investigation 



The presence of complex potentials in these systems makes them ideal candidates to be 
probed using the quantum Hamilton-Jacobi (QHJ) formalism, since this approach has been 
formulated in the complex domain jisl-ll^. In the QHJ formalism, the singularity structure 
of the quantum momentum function (QMF) plays a crucial role in the determination of 
the eigenvalues and the eigenfunctions. Recently, we have studied the structure of the 
QMF in the complex domain for exactly and quasi-exactly solvable Q| potentials. 
In this light, it is extremely interesting to investigate the properties of the QMF of the PT 
symmetric Hamiltonians, in order to systematically find out their differences and similarities 
with exactly solvable (ES) real potentials, as also with the quasi-exactly solvable (QES) ones 
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22j |. It is worth mentioning that, as compared to the solvable potentials, the 



QMF of QES models reveal significant differences in their singularity structure. 

The goal of this letter is to investigate the structure of the QMF of a class of PT symmet- 
ric Hamiltonians, for which the eigenfunctions and the eigenvalues are also simultaneously 
obtained, through the QHJ approach. The differences and similarities of these novel systems, 
with their ES and QES counter parts, are clearly brought out. 



II. QUANTUM HAMILTON-JACOBI FORMALISM 

In this section, we proceed to describe briefly the QHJ formalism and its working. More 
details can be found in our earlier paper With the definition of QMF: 

P = -^^^(InV'), (1) 

the Schrodinger eigenvalue equation Hip = Eip can be cast in the form of the Riccati 
equation, 

p2 _ if^rp' = 2m[E - V{x)]. (2) 
To find the eigenvalues, Leacock and Padgett [3, 16 1 suggested using the following exact 



quantization condition for the bound states of a real potential: 

— </" pdx = nh, (3) 
2n Jc 

where, the contour C encloses the n moving poles in the complex domain, corresponding 
to the nodes of the wave function located in the classical region. This quantization rule 
is an exact one and follows from the oscillation theorem in Strum-Liouville theory. Using 
this quantization rule, Bhalla et. al 0, l^l studied several ES models and showed that the 
eigenvalues could be obtained without obtaining the eigenfunctions. Briefly, the integral in 
Eq.® is evaluated using the knowledge of the singular points of p{x) outside the contour C 
and their corresponding residues. 

The singularities of p consists of fixed and moving ones. The fixed singular points of p 
corresponding to the singular points of the potential, can be found by inspection. From 
Eq.©, one can see that the nodes of the wave function correspond to the moving poles of 
the QMF. In general, the QMF may have other moving poles at locations which are are 
not easy to determine. However, the residues can be computed by substituting a Laurent 
expansion in the Eq.(j2)). In fact, the residue at a moving pole can be easily seen to be 

Application of the QHJ to PT symmetric potentials requires an approach different from 
the one used earlier. In the absence of a generalization of the oscillation theorem j^, for 
this class of potential, it is not clear whether the quantization rule Eq.Q, is valid. Even 
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in the case of the violation of this quantization rule it is not clear which contour should 
be used. Earlier studies indicate that, for all ES and QES models, after a suitable change 
of variables, the quantum momentum function has a finite number of moving poles and the 
point at infinity is at most a pole. In this case the quantization condition still holds for a 
contour, which encloses all the moving poles. We shall assume this to be the case for the 
PT symmetric models to be taken up in this paper. 

When one attempts to compute the residue at a pole using the QHJ, one gets two answers. 
A boundary condition, in the limit h 0, has been suggested by Leacock and Padgett to 
select the right residue. Several other conditions, such as square integrability, have also been 



utilized in earlier papers for this purpose [17|,ll8[. In some cases j27|, 28|, in the absence of 
any criterion to select a residue, one must consider all the values which are consistent with 
the other equations of the theory. 

In studying the complex potentials, we will first perform a suitable change of variable and 
bring the resultant equation to the QHJ form so that the potential is replaced by a rational 
function. We analyze the potential introduced by Khare and Mandal in section 3 and the 
complex Scarf potential V{x) = —A sech — iB tanh x sech x with A > in section 4. After 
comparing the structure of QMF for PT symmetric, ES and QES systems, we conclude in 
the final section, with the remarks about problems and future directions of work. 



III. KHARE-MANDAL MODEL 

The potential expression, for the Khare-Mandal model is given by, V{x) = — (Ccosh2a; — 
iMY- This potential has complex or real eigenvalues depending on whether M is odd or even 
0,1^. It is worth noting that, parity operation in this case is given by a; — > z7r/2 — x, whereas 
the time reversal remains same as the conventional i ^ —i. Using the QHJ formalism, we 
first obtain the QES condition, for the odd and even values of M. Subsequently, the exphcit 
expressions for the eigenvalues and eigenfunctions for the cases, M = 3 and M = 2 are also 
obtained. The QHJ equation in terms of g = ip, setting U = 2m = 1, can be written as 

q'^ + ^ + E+ (Ccosh 2x - iM? = 0. (4) 

dx 

This above form has the advantage that, the residue at each moving pole is one. To bring 
the potential to the rational form we do a change of variable t = cosh2x. Substitution of 
this in Eq.(jll) gives 

■ dq 



q^ + 2Vt^^^ + E + {(t-tMf = 0. (5) 
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One observes that the coefficient of ^ is not one. Hence, in order to bring the above equation 
to the form of Eq. (0]) , we define 



g = 2(Vt^)0, = x--^-^, (6) 



which transforms (0) to 



^ dt ^4(t2-l)2 4(t2_l) ^ ' 

which has the convenient form wherein the residue at a moving pole is one. From here on, 
X will be called as the QMF and Eq.Q, the QHJ equation. 

Singularity structure of x- As already explained in the previous section, we shall assume 
that X has a finite number of moving poles in the complex t plane and that the point at 
infinity is at most a pole. Besides the moving poles, x has fixed poles at t = ±1. It is seen 
from Eq.((7j) that the function x is bounded at t = oo. Assuming that x has only these above 
mentioned singularities, we separate the singular part of x and write it in the following form: 

X-r^H-^-^f^C. (8) 

Here hi and h'^ are the residues at fixed poles t = ±1; P„ is a polynomial of degree n and 
equals 11^=1 — tk), where t^s are the locations of the moving poles of the QMF. C gives 
the analytic part of x and is a constant due to the Liouville's theorem. From Eq.(7), one 
can see that for large t, x goes as which are the values of C. 

To find the residues at the fixed poles, say t = 1, one needs to expand x iii a Laurent 
series: 

X = j^ + ao + ai{t-l) + --- . (9) 

Substituting this in Eq.© and comparing coefficients of different powers of t, one obtains 
the following two values for bi: 

6, = j,i (10) 

Similarly the two values of residues at t = — 1 turn out to be. 

Behaviour at infinity : It has been assumed that the point at infinity is an isolated 
singularity. In order to find leading behaviour of x at infinity, one expands x ^ls, 

X = ao + ^ + ^ + ---. (12) 
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Substitution of Eq. (fT^ which in Eq.Q, gives A = along with, = iy, which is equal 
to C . Due to this A takes the following two values: 



bi+b'+n 
t ' 



This should match with the leading behaviour of x coming from Eq.®, which is 
for large t. Hence equating the two equations, one obtains 

bi + b[+n = \. (14) 

From Eq. pOj) and (fTT|) . we see that the right hand side of Eq. p4|l is positive. Hence for 
Eq.(14) to be true, we choose only the positive value of A i.e, A = which means we 
choose ao = C = It should be noted that, there is no way of choosing a particular 

value of residue at a fixed pole, since one does not have information regarding the square 
integrability of the solutions. Hence, one needs to consider both the values of bi and b[. 
Thus taking all possible combinations of bi and b'l in Eq. (fT!?|l . one obtains the QES condition 
for each combination along with a constraint on M, as given in table 1. From table 1, we 
see that sets 1 and 2 are valid only when M is odd and sets 3 and 4 are valid only when M 
is even. 

Forms of the wavefunction : From Eq.(^ one obtains ip{x) is terms of p with h = 2m = 1, 

as 

^(x) = exp{i J pdx). (15) 
Doing the change of variable and writing p in terms of one gets, 

m =exp/f-^ + -^ + :^ + ^- --^] dt. (16) 

Hence, one can substitute sets 1 and 2 in Eq. ()l(i|l if M is odd and sets 3 and 4 if M is even 
to obtain the form of the wavefunction. The expression for the wave function is in terms of 
the unknown polynomial P„, where n gives the number of zeros of P„. In order to calculate 
the polynomial, we substitute x from Eq.® in ((7j), to get 

6i &'i iC\ bj-bi ib')l-b[ e + 2 



(17) 




t-1 t + i 2 J {t-iy 

E+{Ct- iMf + 86i6'i - AC{e - 1 
4(t2 - 1) 




This leads to n linear homogeneous equations, for the coefficients of different powers of t in 
Pn- The energy eigenvalues are obtained by setting the corresponding determinant equal to 
zero. The explicit eigenvalues and eigenfunctions are obtained for M = 3 and M = 2 cases. 
Case 1: M =3, 
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Here M is odd so we can use sets 1 and 2 from table 1 and get the required results as below. 
Set 1 : bi = i, b[ = i and = 1. 

This implies P„ is a first degree polynomial say Bt + D. Substituting these values in Eq. p7|) 
and comparing various powers of t, one obtains a 2 x 2 matrix for B and D as follows 



1 + -< \ ( B 



4 

a2 



0. (18) 



Equating the determinant of this matrix to zero, one obtains the two values for energy and 
the polynomials as 

E = 7-('±2^1-4C, Pi = |(2t-^(l±V'l-4e)). (19) 

Substituting the values of bi, b[ and Pi in Eq. lfTBj) gives, the two eigenfunctions corresponding 
to the two eigenvalues : 



^jj(x) = eT=°«h2x |'2cosh2a; - ^(1 ± ^1 - 4(2) j . (20) 
set 2: bi = |, b[ = | and n = 0. 

Here we see that n = implies P is a constant. Substituting these values in Eq.(fT7| and 
proceeding in the same manner as before one obtains, 

E = 5 - C^ i^{x) = e^'^^'^^^sinh 2x (21) 



which are the known results LD, |9|| . Below we elaborate on the case when M is even. 
Case 2: M = 2, 

In this case, one makes use of sets 3 and 4 in table 1 and proceed in the same way as was 
done for case 1. 

set 3 : bi = \, b[ = | and n = 0. 

the eigenvalues and eigenfunctions obtained are 

E = 3-C^ + 2zC, ^{x) = e-^''''^^%cosh2x + lY/^. (22) 

set 4 ■' bi = |, b'l = I and n = 0. 
In this case, one obtains 

E = 3-C^ -2iC, ^(x) =e^™^^2"(cosh2x- 1)1/1 (23) 

These match with the solutions given in P| . Thus for any given positive value of M, odd or 
even, one can obtain the eigenvalues and eigenfunctions for the Khare-Mandal potential. In 
the next section, we study the complex Scarf -II potential. 
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IV. COMPLEX SCARF-II POTENTIAL 

The expression for the complex Scarf-II potential is given by 

V = Asech^x + iB sechxtanhx. (24) 

Note that, unlike the previous case, here parity operation is given by a; — > —x, time reversal 
operation remaining the same. The corresponding QHJ equation, in terms of q, where 

-\ V E — Asech "^x — iB sech x tanh x = 0. (25) 

dx 

Carrying out the change of variable, y = zsinhx, proceeding in the same manner as before, 
one obtains the QHJ equation for x • 

2 dx 2 + E A-By _ 

^ ^d^^ 4(1-1/2)2 - " (1-1/2)2 - ' (26) 

where 

X=(0-^^^], q = ^i^J^^^. (27) 

Along with the n moving poles with residue one, x has poles at y = ±1. We assume that 
except for these poles there are no other singularities and that the point at infinity is an 
isolated singularity. The residue at y = 1 and —1 are respectively given by. 



and 



b[ = l±y\ + A + B. (29) 
As in Eq. (12), considering the behaviour of x at infinity one gets, 

A = ^±v^. (30) 

It should be noted that, unlike the previous QES case, energy explicitly enters in A. As 
will be soon seen, this is the reason all the energy values can be obtained here, making this 
exactly solvable model. 

As seen in the earlier section one can write x terms of its analytic and singular parts 

as 

y-l y+l Pn 

where C is the analytic part of x- C* is a constant due to Liouville's theorem, which turns 
out to be zero. The leading behaviour of x fo^' large y from Eq.(|31|) is of the form, 
This coefficient of -, should match with A in Eq.(30) i.e., 



h^ + h\+n = A. (32) 
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This gives the energy eigenvalues as 



-E=ih + b[+n-^y. (33) 
The wave function in terms of x is 

exp(/..(-A_^_^, I (34) 

which is equal to 

^{y) = {y-irP{y + l)-'^Pniy), (35) 

where, —p = 6i — | and —(l = b[ — ^. To obtain the polynomial, one needs to substitute 
Eg . pip in (j^Bj) . which yields the following differential equation 

Pn + + ^) + G(y)P„ = 0, (36) 

where 

(Ajbl -h + bf - b[)) + 1 + AE)y^ + 2yiAibl - &i - bf + b[) + 2B)) , 

Giy) = ^^3Tp + 



(4(6? -bi + bf -b[) + 2-AA- AE) 



(37) 



Substituting the expression for E from Eq. ljH^ in (jH^ . one obtains 

(1 - y^)P" + P'(2(6i - b[) - 2{bi + b[)y) + n(n + 2{bi + b[ - 1) + 1) = , (38) 

which is in the form of the Jacobi differential equation and hence the polynomial Pn{y) = 
p2fei 1,26^ i^^-j Jacobi polynomial. Thus the complete expression for the wave function 

can be written as, 

i:{x) = (isinhx - lf^~^ismhx+lf^'-^P^^^^'^'^^'^'\ismhx), (39) 

which matches with the answer in if bi and b[ are written in terms of p and q respectively. 
Note that in this whole process, we had written the expression of the eigenvalues and the 
eigenfunctions in terms of bi and b'l which have two values. No particular value has been 
chosen. Hence, we need to choose one value of each residue to remove this ambiguity. For 
this case, the solutions are known to satisfy the property ?/'(±oo) —>■ 0. We make use of 
this condition to choose the right values of the residues. For this purpose, we consider 
two conditions on the potential parameter A and B. For each condition, we see that the 
particular values of residues, which are chosen using the property of square integrability give 
physically acceptable results. 
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Case 1: \B\ > A + I 

With this restriction on A and B, the residues at y = ±1 becomes, bi = ^±^JB — A — j 



and 6'^ = I ± A + B + j. In the hmit y — * oo, the wave function in Eq.(35) goes as 

^{y) ^ y^^+b'i'l+^. (40) 

The above equation, with the values of bi and b[ substituted becomes 



For il){y) to go to zero for large y, 6'^ = | + \\J A + i? + ^ is ruled out. Hence the choice of 
the residues for this case will be 



b^ = -±-JB-A-- b\ = ---JA + B + -, (42) 
2 2V 4 ^ 2 2V 4' ^ ' 



with the restriction on n as 



n<]-\ A + B + --]-. (43) 



For these values of the residues at the fixed poles, ip takes the form, 

ip = (zsinh X — l)3^^(isinhx + l)^~^P^*''''*(isinhx), (44) 



where r = J B — A — \ and s = \ A + B + \. The expression for energy is 
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E = -\n + \-\\\A + B + -±iJB-A--\ \ , (45) 
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with the condition on n as 



-<\{i + ^-B + \'Ia + B^\-\. (46) 
Case 2 : \B\<A + \ 

Proceeding in the same way as above, one obtains the choice of the residues as 



6i = - - -\- + A-B, b\=-- -\- + A + B. (47) 

^ 2 2V4 ' ^ 2 2V4 ^ ' 

The wave function is given by, 

iIj = (zsinhx — l)2~2(isinha; + 1) 2^2p^^''^(isinh x), (48) 



where n = Jj + A — B and u = Jj + A + B with the energy 



E = -{n + ---\j\ + A-B + J\ + A + B\ \ . (49) 



Thus the answers match with those given in 
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TABLE I: This table gives the QES condition and the number of moving poles of x for each 
combination of bi and b'l for the Khare-Mandal model 

set bi di n = X — bi — b'l Condition on M QES condition 



^44 2 

9 3 3 M 
^44 2 

o 3 1 M 
"^4 4 2 

4 1 3 M 
^44 2 



1 

2 


M = 


- odd, 


M > 1 


M = 


2n + l 


3 
2 


M = 


- odd, 


M > 3 


M = 


2n + 3 


1 


M = 


even. 


, M > 2 


M = 


2n + 1 


1 


M = 


even. 


, M > 2 


M = 


2n+ 1 



V. CONCLUSIONS 



In conclusion, PT symmetric potentials belonging to the QES and ES class have been 
investigated through the QHJ formalism. The QES solvable Khare-Mandal potential has 
complex or real eigenvalues, depending on whether the potential parameter M is odd or even. 
The singularity structure of the QMF for these two cases is different. For the case when 
M is odd, one observes from table I that the solutions fall into two groups, which consist 
of solutions coming from sets 1 and 2. For a solution belonging to a particular group, the 
number of singularities of the QMF are fixed and consist of both real and complex locations. 
This kind of singularity structure of the QMF has been observed in the study of periodic 
potentials Though the solutions, for M even, fall into two groups coming from sets 3 
and 4, they all have same number of singularities, which again can consist of complex and 
real poles. This singularity structure is same as observed in the ordinary QES models [isf . 

Coming to the case of exactly solvable PT symmetric potential, the location of the moving 
poles can be either real or complex. In the specific example of complex Scarf potential, it 
turns out that all the moving poles are off the real line. In contrast for the ordinary ES 
models the moving poles are always real. For both the cases, the number of moving poles 
of the QMF, characterize the energy eigenvalues. 
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